Direct and alignment-insensitive measurement of cantilever curvature 
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Silicon-based microfabrication has enabled not 
only the electronics revolution but also made 
micromechanics nearly as ubiquitous, with ap- 
plications from motion sensing to biochemical 
analysis [1, 2], The atomic force microscope [3, 4], 
where the motion of a small tip at the end of a 
cantilever traces nanoscale features on surfaces, 
provides the fundamental paradigms for nanome- 
chanical metrology, including optical readouts of 
cantilever displacement and bending. Such read- 
out typically requires careful and costly align- 
ment of mechanical and optical elements. Here 
we demonstrate an alternative approach based on 
the recognition that in many applications we are 
interested not so much in the motion of a tip as in 
the overall curvature of a cantilever. In particu- 
lar, near field imaging of entire cantilevers yields 
diffraction patterns providing precise measures 
of the tilt, curvature and higher order bending 
components of cantilevers. Even while we have 
used very inexpensive components and no care- 
ful alignment is required, we obtain sensitivity to 
nm-scale motion of the cantilever end. 

There are many methods to measure the displace- 
ment of microstructures such as AFM cantilevers [4]. 
They are based on various physical principles, including 
optics [5], piezoresistance[6, 7], field-effect transistors [8] 
and capacitance [9]. The method most commonly used, 
and still the most sensitive and reliable, is the optical 
lever or optical beam deflection technique (OBDT)[10, 
11] -implemented in the AFM market- and optical in- 
terferometry [12-14]. No detection method is optimal 
for all types of measurements and the growing use of 
cantilevers as multimarkers for biomolecules imposes its 
own challenges [15], particularly in the case of arrays of 
several cantilevers. Cantilever arrays are used to obtain 
simultaneous detection of different targets, increase sta- 
tistical significance, in-situ control and for differential 
measurements [16, 17]. They require detection systems 
with a complexity that typically scales with the number 
of cantilevers in the array. Examples of such systems 
include arrays of illuminating lasers with a single multi- 
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FIG. 1. Concepts for standard and proposed optical detection 
methods, a) The optical beam deflection technique (OBDT) 
requires a focused beam carefully aligned to a cantilever and 
the center segmented photo-diode, b) In the No-Alignment 
Nearfield Optical Bending Estimation (NANOBE) technique, 
a broad beam illuminates the whole cantilever array gener- 
ating a diffraction pattern projected onto a CMOS detector. 
The intensity profile of the diffraction pattern is insensitive 
to misalignments but sensitive to the details of the cantilever 
bending. 



plexed detector [18] or 2D scanners of a single laser beam 
employing voice-coil actuators [19], which are simply ex- 
tensions of OBDT for cantilever arrays. 

The main drawback of the common optical techniques 
is the difficulty of accurately aligning each illuminating 
source with its corresponding cantilever and detector, 
making it difficult to measure large number of cantilevers 
without an elaborate pre-measurement protocol. A fur- 
ther limitation particularly of OBDT (Fig la), comes 
from the fact that the observed quantity is the local 
change in angle and displacement of the lever at the 
point of illumination, and therefore cantilever tilt and 
bending cannot be distinguished in a single measure- 
ment. Also, a bending model must be assumed in or- 
der to estimate the true beam curvature. Efforts to 
measure the whole cantilever profile have partially ad- 
dressed the latter issue using scanning or an array of 
Light-Emitting Diodes (LED) [20]. Research has revealed 
that commonly assumed bending profiles may not be as 
realistic as expected[21], providing further justification 
for the development of a simple technique for imaging 
cantilever bending. 

In this Letter we describe a new method for the di- 
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rect and calibrated measurement of cantilever bending 
profiles that requires no scanning or alignment of the 
illumination and does not assume particular bending 
models. We show that an out-of-focus image of the 
whole cantilever array is enough to estimate at least the 
first four coefficients of a polynomial describing the can- 
tilever shape for each of the cantilevers in an array in- 
dependently, making it suitable for multi-marker biolog- 
ical essays. We name this technique NANOBE for No- 
Alignment Nearfield Optical Bending Estimation. 

To detect bending profiles we illuminate a cantilever 
array with a homogeneous monochromatic plane wave 
and measure the intensity pattern of the reflected light 
(Fig lb). Based on the Huygens-Fresnel principle, we 
model the reflected component from the finite-size can- 
tilever as a rectangular source in the plane (£,77). We 
express the wave amplitude in the observing plane (x, y) 
as the convolution integral [22] 
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where [/(£, 77) is the function defining the amplitude 
and phase at the source. We model the shape of the can- 
tilever as a rectangle of dimensions (w,l) and the bend- 
ing profile by a polynomial function P(£) = XI c iC with 
i > 1 ; although here we assume bending only along the 
£ direction which is the long axis of the cantilever, our 
results can be generalized to P(£, 77) taking account of ar- 
bitrary curvature. The bending of the surface introduces 
a difference in the optical path, which introduces a phase 
(j) = 47rA _1 P(£). The field amplitude at the source is well 
described by 
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where rect(x) = 1 for \x\ < 1/2 and zero otherwise. 
We concentrate attention on the x coordinate along the 
cantilever and rewrite equation 1 to separate variables 
U(x,y) = —iex.-p(ikz)X(x)X(y) and considering that the 
observed intensity I(x,y) = \U(x,y)\ 2 = \T(x)\ 2 \I(y)\ 2 
we analyze the value of 
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FIG. 2. Intensity map I(x) for a cantilever as a function of 
different figures of interest. Insets are vertical cross-sections 
of the intensity map showing the intensity profiles for specific 
abscissa values. The horizontal axes are: a) Fresnel number 
F — L 2 (Az) _1 , where L is the cantilever length, and z the dis- 
tance from the cantilever at which the intensity is measured, 
b) c\ the cantilever tilt, c) C2 the cantilever curvature and d) 
C3 the cubic bending. 



with m = 1 + 4c2Z , s = 2c\z and 



exp 



[ik (2 



I c 2 x +£j_C 



(5) 



If we neglect cubic and higher order terms in the poly- 
nomial description of the cantilever curvature i. e. assum- 
ing Ci « Mi > 2, we can introduce a change of variables 
and a corresponding change in integration limits to find 
a familiar result. Using 
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We now use the first two terms of P(£) to complete a 
squared binomial for £ and for convenience rewrite equa- 
tion 3 as 
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we obtain the solution for a rectangular slit 
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where C(a) and S(a) are Fresnel Integrals de- 
fined as C(di) = Jq Z cos (7ra 2 /2) da, and S(oti) = 
(4) J** sin (ira 2 /2) da and the integration limits a\ = 
a(-L/2) and a 2 = a(L/2) 
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The inset of figure lb shows calculated 2D patterns given 
by I(x, y) = \X(x)\ 2 \X(y)\ 2 resembling experimentally 
observed patterns. Figure 2a shows I(x) for a rectangu- 
lar slit or a flat cantilever i.e. q = Mi as a function 
of the Fresnel number F = L 2 (Az) _1 , where L is the 
cantilever length, and z the distance from the cantilever 
at which the intensity is measured (insets are intensity 
profiles for specific abscissa value). In what follows we 
analyze the influence of each coefficient q / on the 
observed diffraction pattern. 

Form of equation 6 itself and the figure 2b we see that 
a non vanishing tilt c\ results in a shift s = 2c\z of the 
diffraction pattern, as expected from the tilt of any re- 
flecting surface. Similarly from figure 2c, provided z ^ 
a non-zero c 2 causes a magnification m = 1 + 4c 2 z in the 
pattern size and m _1 in the intensity. The magnification 
does not come from the negligible displacement of the 
cantilever ends, but is intrinsic to the curvature of the 
reflecting surface. Both effects do not otherwise distort 
the shape of the diffraction pattern. 

Higher order contributions to the curvature are more 
easily studied by numerically integrating equation 4. Fig- 
ure 2d shows the Fresnel diffraction pattern for cantilever 
deflected as P(x) oc x 3 . Cubic bending shown in figure 2d 
feature intensity gradients and pattern shifts, both pro- 
portional to C3. Quartic bending (not shown) creates a 
curvature in the intensity profile as well as a change in 
the pattern length. 



I. EXPERIMENTAL 

Figure 3a illustrates the experiment where a broad 
(16.8mm) collimated laser beam (A = 660nm) illumi- 
nates a cantilever array. A 4X microscope objective and 
a CMOS sensor focused in a plane at a distance z from 
the cantilever maps the diffraction pattern from the can- 
tilevers. The array (purchased from IBM) consisted of 
up to eight rectangular silicon (100) cantilevers, each 
measuring 500/im in length, 100/im in width, and 0.9/im 
in thickness and a nominal spring constant of 0.027V/m. 
The cantilevers were coated on one side with a 2nm tita- 
nium adhesion layer and then 20nm gold. The diffraction 
pattern size is calculated in units of pixels by software 
simply by using the distance between intensity threshold 
crossings or, more preferable, using the second central 
moment of a intensity profile (See Supplementary 
Information) . 
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FIG. 3. Experimental measurements using NANOBE 
method, a) Setup schematics: A cantilever array is illumi- 
nated with a broad laser beam through a cubic beam-splitter 
and the reflected light captured by a lens and CCD. b) A sin- 
gle cantilever is pushed with a glass tip mounted over a cali- 
brated piezoelectric actuator, c) Diffraction image (640 x 528 
px) of the pushed cantilever (top pattern) are shorter than the 
pattern from the relaxed cantilever by a factor of (1 + 4bz) 
d) The experimentally measured pattern size changes linearly 
with the cantilever tip displacement, e) The curvatures of 
six cantilevers as a function of time. Cantilevers functional- 
ized with MHA (top) and HDT (bottom) have different static 
bending and different sensitivity to pH changes. 



We test the method experimentally for deflections in 
excess of one micron by measuring the diffraction pattern 
as a function of cantilever bending where we poke the 
cantilever end in air with a glass tip (Fig 3b) attached to 
a well-calibrated piezoelectric actuator (P-363 PicoCube 
XYZ Piezo Scanner by Physik Instrumente GmbH). The 
concentrated load at the free end is expected to cause a 
bending profile proportional to x 2 (3L — x), where L is 
the cantilever length [23]. Along with the controlled dis- 
placement we observe both a change of the size of the 
diffraction pattern as well as a gradient in the intensity, 
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FIG. 4. The cantilever bending caused by small temperature 
changes are easily resolvable. The temperature protocol was 
set to 25 o C,25.05 o C,25.10 o C,25.15°C and 25.20°C. The light 
trace corresponds to the raw data and the darker trace is a 
10 point moving average. 



characteristic of both parabolic and cubic bending (Fig 3c 
top pattern). The change in pattern size displays the pre- 
dicted linear relationship with the displacement caused 
by the actuator (Fig 3d). 

We further test the method in fluid by functionalizing 
the surface of individual cantilevers with either mercap- 
tohexadecanoic acid HS(CH2)i5COOH, here abbreviated 
MHA, or hexadecanethiol HS(CH2)i5CH3, here abbre- 
viated HDT via incubation in an array of glass micro- 
capillaries. The protonation/deprotonation of MHA at 
different pH causes a differential surface stress inducing 
cantilever bending. HDT has identical chain length and 
similar packing density to MHA but differs in the ter- 
minal methyl group which is non-ionizable and therefore 
is used as reference [24]. The cantilever array is then 
mounted in a liquid cell with a sapphire window to al- 
low illumination and imaging. An automated system of 
syringes and valves was developed to control the deliv- 
ery of sodium phosphate mono and dibasic solution pH 
4.8 and pH 9.0 at constant rate of 43/iL/min and con- 
trolled temperature of 25.00 ± 0.01°C. The supplemen- 
tary materials show a schematic as well as photograph of 
the instrumental prototype which we have constructed 
for these measurements. Figure 3e shows the distinctive 
changes in pattern size as we alternately flow solutions 
with pH 4.8 and pH 9.0. The main changes in curva- 
ture are attributed to the differential stress caused by 
protonation/deprotonation. The mean relative deflec- 
tion for MHA cantilevers is calculated to be 249. 8nm, 
consistent with previous results [25]. The mean relative 
deflection for HDT is somehow smaller than previous re- 
sults (137. 2nm) giving a differential stress between MHA 
and HDT-coated cantilevers of around 0.0237V/m, almost 
two times bigger than expected [24] . Observing the initial 
curvatures in figure 3e, we see that NANOBE reveals that 
HDT and MHA functionalization causes different static 
bending of the cantilevers as well as different sensitivi- 
ties to pH. The variation in static and dynamic bending 



among cantilevers with the same functionalizations re- 
veals the difficulties in preparing consistent surface mod- 
ification in small areas. The inconsistency upon func- 
tionalization of nominally identical cantilevers remains 
an open question to the authors and is beyond the scope 
of this work, but it should be noted that this information 
is typically unavailable on OBDT measurements whereas 
evident via NANOBE. 

To demonstrate the resolution of small deflections in 
liquid we step the temperature of the cantilever in in- 
crements of 0.05 ± 0.01 o C around 25°C. The different 
thermal expansion of silicon and gold causes the micro- 
structure to bend. Figure 4 shows that the measured 
deflection is well resolved even for these small changes in 
temperature. The inset on the right hand side shows the 
density of measured deflections, featuring well resolved 
maxima and negligible overlap of the individual distribu- 
tions. The cantilever is found to bend around 150nra/°C, 
in good agreement with previous measurements and mod- 
els [26]. 

II. DISCUSSION 

The main output of the NANOBE method is the size 
of a diffraction pattern which is independent of the region 
of the detector that is used to capture it, and therefore 
insensitive to small misalignments of the light source or 
detector and independent of the orientation of the can- 
tilever. The pattern size measured in units of pixels is 
directly translated to an absolute value for the curvature 
using 

c 2 = (So - S c )/(4zS c ) (11) 

where S Q is the observed pattern size in pixels, S c = 
pL Cl where L c is the cantilever length and p resolution 
of the sensor in pixels per length, and z is the distance 
to focus. In contrast to the standard laser beam de- 
flection method where the optical lever arm needs to be 
carefully fixed and measured, our method is entirely self- 
calibrated, i.e. the known geometry of the cantilever ar- 
ray gives the conversion ratio between length and pixels, 
and the cross section of the diffraction pattern reveals 
the effective value of distance z. Nonetheless, if there 
are any doubts, these values are also available to the ex- 
perimentalist by physically measuring the geometry of 
the experimental setup. The changes of the tilt c\ are 
also available from the centroid of the diffraction pat- 
tern, and the higher orders of curvature cs and above, 
can be estimated from the overall intensity gradients for 
each pattern. 

Provided the Fresnel number is large enough, the 
diffraction patterns will be localized with negligible cross 
talking between adjacent cantilevers, therefore allowing 
independent analysis of simultaneously imaged patterns. 

The resolving power of the system is proportional to 
the magnitude of z, to the number of pixels covered by 
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the pattern and inversely to the noise level of the de- 
tector. Our experiments were performed with a sim- 
ple USB 1280 x 1024 px Monochrome CMOS (Thorlabs 
DCC1545M) where each pattern covered only around 400 
pixels, therefore a much improved resolution can be ex- 
pected from a high-end detection system. 

III. CONCLUSIONS 

We have invented a novel near-field method for mea- 
suring cantilever bending which is much more robust and 
cheaper to implement than the standard optical beam 
deflection methods. It relies on curved mirror diffraction 
and is insensitive to misalignment, opening the way to a 
variety of inexpensive - from both manufacturing and op- 



erational points of view - devices for nanometrology. The 
devices range from force microscopes to biochemical as- 
says where the cantilevers could even be weakly tethered 
(or untethered) to substrates. 
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